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ABSTRACT 

We present observational results of a type II burst associated with a CME - CME interaction observed 
in the radio and white-light wavelength range. We applied radio direction-finding techniques to 
observations from the STEREO and Wind spacecraft, the results of which were interpreted using 
white-light coronagraphic measurements for context. The results of the multiple radio-direction finding 
techniques applied were found to be consistent both with each other and with those derived from the 
white-light observations of coronal mass ejections (CMEs).The results suggest that the Type II burst 
radio emission is causally related to the CMEs interaction. 

Keywords: solar-terrestrial relations - Sun: coronal mass ejections (CMEs) - Sun: radio radiation. 



1. introduction 

Over periods of increased solar activity, several coro- 
nal mass ejections (CMEs) can be launched by the same 
or nearby active regions ( Gopalswamy et al.|2QQ5( ) . Dur- 
ing these times of high activity, one or more of these 
CMEs may interact while propagating through the in- 
terplanetary medium. Almost a decade after th e rst ob 



servations of CME - associated shock regions (Burlaga 
^t al. 1987), CME-CME interactions were observed, at 
long wavelen gths and in whi t e-li^ht coron agrap hic im- 



ages, by [Gopalswamy et al. (2001 2002) and Gopal 



swamy 
by the 



(|2004p . 'I'he radio observations were obtained 

Radio and Pla sma Wave Experiment (WAVES, 

[Bougeret et al. 11995 ) on board the Wind spacecraft, 
while the white-light observations were obtained by the 
Large Angle and Spectroscopic Coronagraph (Brueckner] 
|et al^l995) on board the Solar and Heliospheric Observa- 
tory mission. Based on the observational characteristics 
of the CM Es from white- light coronagraph and radio ob- 
servations, Gopalswamyl ([2004 ) concluded that the type 
II radio emission is enhanced and modified due to the in- 



teraction between two CMEs. Gopalswamy et al. (2001) 
suggested that the observed radio enhancements result 
from the increased density in the upstream medium that 
reduces the Alfven speed, thereby increasing the Mach 
number of the shock. This is in agreement with results 
from numerical simulations, confirming that the radio en- 
hancement was lik ely to be pr oduced at the in t eraction 
region shock (e.g. 



Vandas and Odstrcil 2004). Gopal- 



swamy et al. (2001) also mentioned additional possibil- 



ities tor electron acceleration, such as reconnecti on be- 
tween the two CMEs (see also | Gopalswamy | |2004 ) . 
White-light (WL) imagers such as those on board 
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the Solar TEr restrial RElations Observatory {STEREO; 



Kaiser et al. 2008) allow us to observe CMEs out to 



1 AU with the ISun - Earth Connection Coronal and He- 



liospheric Investigation (SECCHI; [Howard et al. 2008) 



The Corl and Cor2 coronagraphs, along with the helio- 
spheric imagers (HIl and HI2) on board, observe CMEs 
at visible wavelengths as they propagate through the he- 
liosphere. Coronograph observations typically reveal the 
three part configuration of CMEs: the bright, dense core 
that is thought to be the erupting filament; the dark, low 
density cavity surrounding the core; and the bright front, 
or the leading edge. 

The stereoscopic observations of the STEREO mission 
allow us to determine the locati on of different CME fea- 
tures in three dimensions (e.g., Frazin et al. 1 2 009 As- 



chwanden and Wulser|20rT||Liu et al. 



et al.||2010D. A similar approach can be used to de- 
termme the position of radio sources in what is called 
radio direction-finding. Several direction-finding tech- 
niques have been implemented using observati ons made 
either by spinning spacec raft like Wind (e.g., Fainberg 



et al. 1972 Reiner et 



jt all 
JTER 



1998) or three-axis s tabilized 
r K) or Cassini (e.g., 'Cecconi 



spacecraft s uch as S'l 

|et al.||2008| [S antolik et" aL ]|2003D. The stereoscopic ca- 
paBility of gri^'A^A'6»/ WAVE5"i |BQu'geret et al.|2008D can 

be used to triangulate the three dimensional position of 
a radio source at a particular frequency, provided both 
spacecraft observe the same source quasi-simultaneously. 
This process can be repeated for different frequencies. 
This technique has been applied with great succes s in 
the past in the study of t ype HI emission ( Gurnett e t al.| 
1978 Reiner et al.||2009 ) but rarely in the study of type 



11 bursts. Several successful campaigns were undertaken 
using a combination of spacecraft, such d,^ Helios, Ulysses 
and Wind (e.g., Baumback et al. 1976 [Reiner et al.| 
1995), demonstrating the success oi direction- finding, by 



mapping the path of accelerated electrons during type 
HI bursts. 

In this paper we study the relationship between the 
interaction of two CMEs and the location of the as- 
sociated radio sources during an event which occurred 
on 2010 August 1. We made use of three space-based 
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instruments with direct ion- findi ng capabilities , namely 
STEREO /WAVES experiments ( iKaiser et al.|2 QQ8; Bale] 
|et al.||2QQ8|) and Wind^ along with white- light data from 
the SECCHI suite onboard STEREO. 

2. OBSERVATIONS AND ANALYSIS 

The period 2010 July 31 to August 2 was character- 
ized by increased solar activity, exhibiting small flares, 
filament eruptions an d coronal mass e j ectio ns (Schrijver 
[and Title||2Qll[ [Temmer et aL||2011[ |Liu et al. 201 1|). 
Of particular interest here is the time during which two 
CMEs (one slow(CMEl), erupted at 02:00UT and one 
fast(CME2), erupted at 07:00UT) interacted with each 
other, resulting in a low frequency type II radio burst 
observed on 2010 August 1 at about 09:00 UT). The two 
CMEs in question can be seen in Figure [l] (left panels 
- Cor2 A, right panels - Cor2 B). Each row corresponds 
to a time either before (top and middle rows) or during 
(bottom rows) the interaction period. 



eclipti c plane of STEREO-A and -5, as described in |Liu| 
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Figure 1. White-light observations of the two CMEs labeled as 
CMEl and CME2 at three different times during the evolution of 
the event. The interaction between the CMEs is clearly seen in the 
bottom frames. 

From these white-light observations, we determined 
the velocity and time of interaction of the expand- 
ing CMEs. An elongation map was constructed from 
running-difference images of Cor2 and HIl along the 



et al. ( 2Q1Q| ). The filamentary structures in the elon- 
gation map (Figure [2]) are the propagating CMEs ob- 
served in the period 2010 August 1-3. Figure [2] shows a 
fast CME (CME2) that intersects and overtakes a slow 
CME launched earlier (CMEl); the region of interaction 
is shown by a dashed box. The average speed of the fast 
CME, derived from Cor2 observations is ^1138 kms~^ 
with a liftoff time of -07:48 UT from the Sun. The liftoff 
time of the slow CME was calculated to be 02:48 UT with 
an average propagatio n veloc ity of 730 kms~^ in Cor2 



fibmmer et aL|2Qll[ [Liu et al.||201l| ). 
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Figure 2. Time -elongation maps constructed from running- 
difference images of Cor2 and HIl along the ecliptic plane for 
STEREO A and B. The dashed box shows the interaction region 
between the two CMEs labeled CMEl and CME2. 

Unfortunately, the STEREO-B Corl and HIl had a 
data gap of about 18 hours starting at 9:20 UT, which 
restricted our analysis. It is important to mention that 
the speed estimates presented above were obtained using 
an algorithm that fits the distance with a linear model. 
These results were then compared to the speed derived 
from adjacent distances with a three-point Lagrangian 
interpolation, obtaining similar results. From the La- 
grangian algorithm, it is possible to derive the error of 
the computation giving ± 315 kms~^ and ± 206 kms~^ 
for the fast and slow CMEs, respectively. It should be 
noted that such large errors are present in all methods 
that measure distances and are not only present in the 
triangulation method used in this study. 

2.1. Radio emission 
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The type II radio burst of interest and its modification 
by CME interaction is shown in Figure [3j This event was 
detected simultaneously by both STEREO and Wind 
spacecraft. The radio emission was characterized by a 
slow drifting feature, first observed by STEREO-B at 
about 9:10 UT and ending around 11:30 UT. The drift 
velocity observed by Mn(i/WAVES, starts at 2000 kHz 
and ends at about 700 kHz. 

As Figure [3] shows, the drifting feature in the radio 
spectra appears to split into two bands at about 9:50 UT. 
This may be related to the properties of the ambient 
plasma or, as we show here, the signatures of two in- 
teracting CMEs. We fit the STEREO /WAVES radio 
spectra to determine the propagation velocity and the 
distance at which the emission could be produced, as- 
sumin g an interplanetary density model (|Leblanc et al.| 
1998). The radio spectrum was divided inTtwx) regions. 



termed upper and lower branches following the observed 
division in the spectrogram. Table [T] shows the radial 
distance to the Sun derived from the model. From the 
model the drift velocity was calculated (Table [T]), reveal- 
ing that the lower branch has a higher drift velocity, indi- 
cating that this part of the type II radio emission was as- 
sociated with the fast CME (CME2). The upper branch 
demonstrated a slower drift velocity which is consistent 
with original velocity of the first and slower CME, indi- 
cating that this was the signature of the slow CME or 
merged CME front. 



Distance (AU) 



Branch 


STEREO-A 


STEREO-B 


Lower 


0.025-0.033 


0.022-0.036 


Upper 


0.025-0.027 


0.025-0.027 


Velocity (kms ^) 


Branch 


STEREO-A 


STEREO-B 


Lower 


1370 


1600 


Upper 


290 


400 



Table 1 

Estimated radial distances to the Sun and drift velocities derived 
from Leblanc et al. ( 1998 ) density model for both branches in the 
radio spectra and for each STEREO spacecraft. 



Since the velocities are determined using an interplan- 
etary density model, the results are highly sensitive to 
any change of the 1 AU electron density (in the model) 
used for this computation. An uncertainty estimate of 
100-200 kms~^ was ob tained by calculating the speed 
using the Leblanc et al.| (p-998) density model for variety 
of ambient electron density values which range from 4 to 
7 cm~^ (as observed in in situ data 4 days after the event 
by Wind). These velocities are comparable to the prop- 
agation velocities derived from Cor2 observations within 
the errors of the measurements, suggesting that the split 
branches observed at 09:50 UT in radio spectra are the 
signatures of the two interacting CMEs. 



2.1.1. Direction-finding 

This event was observed by three spacecraft which, 
with radio direction-finding capabilities, gives a unique 
opportunity to study and locate the region or regions 
responsible for the radio emission in the interplanetary 
medium. There are different techniques that allow us 



to determine the distance at which the emission was pro- 
duced relative to the observer. Some of these make use of 
electron density models, as demonstrated in the previous 
section, which provide a direct correlation between the 
observed freque ncy and the distance (h eight) at which 
they occur (e.g., Leblanc et al.|1998 Reiner et al ^|2007 ) 



However, these techniques do not take into account m- 
homogeneities that may occur in both the interplanetary 
space and/or the ejected material. Also, the propaga- 
tion direction cannot be determined by a density model. 
Other "direction-finding" techniques, which locate the 
region of emission by triangulating the position of the 
radio source at distances of ^0.1-0.5 AU, have been de- 



veloped during the last four decades (e.g., Fainberg et al. 


19721 ICecconi et al. 


20081 Santohk et al.|2003| |Martinez- 


Oliveros et al.||2012 


). 



tion algori thms (M artinez-Oliveros et al.||2012[ Santohk 



et al.|[2QQ3| respectively) to determine the arrival direc 



tion ol radio waves in the frequency range of the High 
Frequency Receiver 1 instrument onboard STEREO. For 
Mnc?/WAVES data, a modulation technique was app lied 
t o re trieve the radio waves' direction of arrival ( Fain berg] 
|et al.||l972| in the range of the Radio Receiver Band 2. 
I'he direction of arrival was then characterized by unitary 
vectors, defined by the azimuths and elevations found by 
the direction-finding procedure for all observations in the 
time range. The spatial positions of the radio sources 
in interplanetary space were fou nd using a geometrical 



triangulation algorithm based on 
simplicity, we will refer to the com' 



Liu et al. 



(2010). For 



Dined direction-finding 
and triangulation techniques as "direction- finding" . 

In order to determine the location of the type II ra- 
dio burst relative to the CMEs, white-light images from 
SECCHI were used. Data from the Cor2 instrument on 
board STEREO B prior to and at the time of the radio 
burst were used to compare the projected radio direction- 
finding results (see Figure [4| . This type of comparison 
has been applied before in the study of type II bursts 
(e.g., Wagner 1982). The results of our analysis suggest 
a close relationship between CME - CME interaction re- 
gion and the type II radio burst. T his agrees with the 
findings of |Gopalswamy et al. (2001 ) who reported a sim- 
ilar result lor a ditterent event. 

The dual views of the STEREO spacecraft were ex- 
ploited to identify the location of the radio burst in 3D 
space, relative to the CME structure. Figure [5] shows the 
front of the CME in heliographic coordinates obtained 
using the Solar Soft package suite (Freeland and HandyJ 
1998). The three dimensional locations ol two CMEs 
are shown in Figure [5] the position of the slow CME 
launched at about 02:00 UT is shown as yellow dots, 
while the position of the fast CME launched at approxi- 
mately 07:00 UT is represented by blue and green dots. 
Two evolutionary times are shown, 09:00 and 10:00 UT, 
as during this period the interaction between the two 
CMEs occurred, with the fast CME overtaking the slow 
one slightly after 09:00 UT (see Figure [2|. Note that 
these three-dimensional observations, shown in Figure [5) 
are projected on the ecliptic plane. 

Figure [s] (left) shows the triangulated position of the 
radio source for the three combinations of STEREO-A^ - 
B and Wind spacecraft at the time of the most prominent 
peaks in the radio flux for three frequencies (925, 975, 
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Figure 3. STEREO- A , STEREO-B and Wind dynamic spectra of the 2010 August 1 type II burst from 09:00 UT to 12:00 UT. The 
plotted frequency range from 125 kHz to 16.025 MHz for STEREO and Wind. The color shading represents the intensity of the radio 
emission measured in arbitrary units. 

and 1025 kH^ see Figure [§. 

The apparent misalignment between the triangulated 
locations from different spacecraft can be explained by 
understanding that the triangulation algorithm searches 
for the position in space where the vectors intersect. This 
intersection does not necessarily occur at the front of the 
emitting region, or its centroid (Figure [sl right frame) . 
Here, it is likely that the source observed is highly ex- 
tended and complex. Therefore, each spacecraft identi- 
fied different regions of the extended source (blue dashed 
ellipse) due to, e.g. the structure of the region and the 
surrounding local plasma density. Another possible ex- 
planation for the apparent discrepancy is that a dense 
region was located somewhere between the type II radio 
source and STEREO- A. This region could scatter the ra- 
dio waves, leading to an apparent shift of the line-of-sight 
source position. 

2.1.2. Time- of -flight analysis 

We examined the timing of the radio profiles at the 
three spacecraft as a control technique to validate the lo- 
cations and results determined by direction-finding. We 
do this by first computing the distance to each spacecraft 
from the extrapolated locations. Then, times of fiight for 
each spacecraft are computed assuming that the radio 
emission travels in a straight line from the source cen- 
troid to the spacecraft at a constant velocity (the speed 

^ For context, the Parker spiral is plotted in Figure [s] and was 
calculated using the formula cf) = (f)o — (GQ/Vsw)r, where r, is 
the radial distance to the Sun, is an arbitrary angle, ^0 is the 
rotational velocity of the Sun (2 kms"-*^), and Vsw is the solar wind 
velocity (400 kms~^). 



of light). The difference between these two times is com- 
pared with the time shift between radio fiux profiles at 
the three spacecraft (see Figure |6| to determine whether 
they are consistent with the source locations found using 
the direction-finding method. This "time-of-fiight anal- 
ysis" assumes that the onsets of the signals at the two 
spacecraft are the signature of radio emission simultane- 
ously emitted from a single compact source. The limita- 
tions of this "time-of-fiight" analysis are in the temporal 
resolution of the measurements and errors inherent in the 
assumptions of compactness and simultaneity. 

We find that the time shift (delay) computed from the 
direction-finding results ranges from ^2 minutes at the 
lowest frequencies decreasing to ^1 minute at the high- 
est, while the observed delay between the peaks of the 
emission received by STEREO- A and -B ranges from ^1 
to ~0 minute, respectively. In our analysis we also make 
use of Wind data. Comparing the times of arrivals at 
the Wind and the STEREO spacecraft computed from 
the direction-finding results, we find an average time de- 
lay between ?^5.2 minutes and ~1 minute, while the ob- 
served time delay between the radio signatures is about 
^2 minutes. We found that the direction-finding and 
the time-of-flight analysis results are consistent within 
the errors inherent to both techniques. The geometrical 
configuration suggested by the time-of-fiight analysis is 
consistent to the one obtained by the direction-finding, 
in which Wind is located closer to the radio source that 
either of the STEREO spacecraft, and also that the radio 
source is located almost at the same distance from each 
STEREO spacecraft. 



A CME-CME Interaction, and its radio and white-light manifestations 



5 



1.5x10^ 



1 .0x10^ 



5.0x10^ 



l-Aug-2010 09:09:04.138 UT 



-5.0x10^ 



.0x10^ 




1 .5x10^ 



1.0x10^ 



5.0x10^ 



-1.0x10^-5.0x10^ 5.0x10^ 1.0x10^ 1.5x10^ 
X (arcsecs) 



l-Aug-2010 10:09:04.153 UT 



-5.0x10^ 



-1.0x10^ 




-1.0x10^-5.0x10^ 5.0x10^ 1.0x10^ 1.5x10^ 
X (arcsecs) 

Figure 4. STEREO-B observation of the 2010 August 1 
CME from Cor2 at two representative times of type II 
radio burst, 09:09 UT and 10:09 UT, with Une-of-sight 
direction-finding results from STEREO-B /WAVES overplot- 
ted in color, where color represents different frequencies. The 
solid white line in lower panel shows the contour of the ex- 
panding CME. 

3. CONCLUSIONS 

During 2010 July 31 -August 2 a series of CMEs and 
their associated type II and type III radio bursts were 
observed. In particular, an interplanetary type II burst 
was detected by instruments onboard the STEREO and 
Wind spacecraft on 2010 August 1 at about 09:00 UT. 
The close timing between the type II radio burst and the 
interaction of two coronal mass ejections suggests that 
the radio emission is a consequence of this interaction. 
A sin i ilar ev ent was analyzed previously by Gopalswamyl 
et al. (2001 ), who concluded that the interaction between 
slow CME and a fast one resulted in the enhancement of 
the radio emission during the transit of the fast CME 
shock front through the core of the slow CME. 

Using white-light and radio observations we estimated 



the propagation velocities of the two CMEs. We 
found that the velocities derived from radio observa- 
tions are comparable to the propagation velocities de- 
rived from coronagraph observations. This suggests that 
the branches in the radio spectra, observed at 09:50 UT, 
are the signatures of the two interac t ing C MEs. Using 



the density model of Leblanc et al. (1998) we also eS' 



timated the distance at which the radio emission was 
produced, was between 0.025 and 0.043 AU. This is in 
agreement with the radio direction-finding results, which 
give a distance about 0.01-0.05 AU. From white-light 
observations, we determined that the shock front prop- 
agated ~20° east of the Sun -Earth line (i.e. between 
STEREO-B and Wind), which is about the same angu- 
lar separation derived by the direct ion- finding technique. 
The obtain ed prop agation direction i s in agreement w ith 
finding of Temmer et al. (2011) and Liu et al. (2011). 
We successfully applied three radio direction-finding 



techniques (Fainberg et al 
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II radio burst and determmed the direction of arrival of 
the radio emission. The data analysis shows that the 
radio sources locations are spread over a large area cov- 
ering about 4°, suggesting that the radio source has an 
extended and complex structure in nature, perhaps com- 
posed of multiple radio emitting regions which may have 
a common origin. We found good consistency between 
the triangulated white-light positions and the Wind- 
STEREO-B triangulated positions. Using STEREO-A, 
we found a discrepancy that can be explained by the com- 
plexity of the source and the surrounding material. Since 
neither the emitting region nor the medium are homoge- 
nous, it is possible that the radio source was partially 
occulted in the direction of STEREO-A by a dense solar 
wind region. This may explain the relatively low power 
observed in the STEREO-A spectrogram and can also ac- 
count for scattering of radio waves, which consequently 
will shift the apparent position of the radio source. 

By comparing these positions with white-light features 
in the STEREO coronagraph d ata an d their derived po- 
sitions as described in Section [2.1.1[ we found that the 
radio emission is the result of the interaction between 
two expanding CMEs. Figure 4] shows that the positions 
derived from the direction-finding match the features ob- 
served in the coronagraph images, suggesting the relation 
between the type II radio emission and the interaction 
region of two expanding CMEs. 

Radio direction-finding has proven to be a powerful 
technique in the study of CMEs and associated type II ra- 
dio bursts. By using these techniques, it is possible to de- 
termine the heliographic distance of a radio source, which 
is independent of any density model. In general, the lim- 
itation of these techniques is given by the frequency of 
observations and the properties of the radio emission re- 
gion. In the case of metric wavelengths this error can 
be about 1° in azimuth and elevation. Nevertheless, our 
study shows that with good observations, the techniques 
give results that can be directly compared with observa- 
tions at other wavelengths and show the likely emission 
region. The application of radio direction-finding meth- 
ods to data acquired by future missions, such as Solar 
Probe Plus and Solar Orbiter, will prove to be crucial in 
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Figure 5. Left: location of the geometrically triangulated positions of the radio sources in interplanetary space for all operational 
frequencies as seen from the top. The Parker spiral is plotted in gray for context. Right: radio emission scenario, showing the 
possible emission region as an extended source propagating between STEREO-B and Wind. The dots represent positions of the 
two associated coronal mass ejections at different times. The red symbols show the intersection between line-of-sight vectors 
from the spacecraft represented by the arrows. These are results projected on the ecliptic plane. The overall regions where the 
direction-nding positions are located are represented by the color shaded areas. The dashed ellipse shows the area covering all 
direction results. 



our understanding of CMEs and type II radio bursts. 
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